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Electric flux tubes are a well known attribute of the QCD vacuum in which they manifest con- 
finement of electric color charges. Recently, experimental results have appeared suggesting that not 
only those objects persist at temperatures T ~ near the QCD phase transitions, but their decay 
is suppressed and the resulting clusters in AuAu collisions are larger than in pp (i.e. in vacuum). 
This correlates well with recent theoretical scenarios that view the QCD matter in the T ~ Tc region 
as a dual-magnetic plasma dominated by color-magnetic monopoles. In this view the flux tubes are 
stabilized by dual- magnetic currents and are described by dual-magnetohydrodynamics (DMHD). 
In this paper we calculate classically the dissipative effects in the flux tube. Such effects are associ- 
ated with rescattering and finite conductivity of the matter. We derive the DMHD solution in the 
presence of dissipation and then estimate the lifetime of the electric fiux tubes. The conclusion of 
this study is that a classical treatment leads to too short of a lifetime for the flux tubes. 

PACS numbers: 



I. INTRODUCTION 
A. Motivation 

Non-perturbative phenomena in both the QCD vac- 
uum and in the finite temperature/density QCD mat- 
ter have been the subject of intense studies for a long 
time. Various phenomenological approaches have been 
proposed during the last few decades to tackle them. One 
simple example is the stochastic QCD vacuum model [!_ 
which provided a picture of the gluonic " condensate" and 
correlations as a constant fields, another is the instanton 
liquid model [5] that explained several phenomena related 
to the SU{Np) and the C/a(1) symmetry breaking. 

Most researchers in the field hold to a generic picture 
of confinement introduced by t'Hooft and Mandelstam 
in the 1980's: the so-called "dual superconductivity" . It 
suggests that confinement is a dual Meissner effect, and 
that a "coil" which prevents the color-electric flux tube 
from spreading out is produced by a magnetic supercur- 
rent. The condensate of Cooper pairs produced in usual 
superconductors is substituted by a Bose-Einstein con- 
densate of some objects possessing color magnetic charge. 
Presently, however, there is no concrete understanding 
of the field configurations that are accountable for color 
confinement. The effort for identifying the correspond- 
ing topological objects responsible for color confinement 
is still continuing to this date. Monopoles, dyons, or 
their composites, induced by fermions, have been shown 
to provide a mechanism for confinement and chiral sym- 
metry breaking in the Af=2 SYM theory P]. But it is 
still unclear whether such objects can be concretely car- 
ried into the QCD-like theories without involving scalars. 
The seminal paper by Polyakov '121 set an early example 
of confinement in (2+l)D Yang-Mills theories by means 
of gluomagnetic monopoles. Recently, Unsal |4] has pro- 
posed an interesting (3-|-l)D extension of the latter model 



of confinement for QCD-like theories, involving fermions, 
by invoking certain composite objects endowed with total 
magnetic and zero topological charges. However, here too 
it remains to be seen whether such objects can be truly 
responsible for confinement, particularly in the frame- 
work of lattice gauge theories. 

In the last few years, this subject has been revital- 
ized by experimental studies at the Relativistic Heavy 
Ion Colhder (RHIC), where hot QCD matter is produced 
and studied. One conclusion stemming from the analysis 
of RHIC results is that the quark-gluon plasma (QGP) in 
the l-2Tc temperature interval is a strongly coupled fiuid 
(sQGP). Liao and Shuryak [S] (a review on the subject 
can also be found in [6|) have related this finding with the 
so called "magnetic scenario" , arguing that while color- 
electric objects (quarks and gluons) interact stronger and 
stronger when approaching the critical temperature Tc, 
the Dirac condition implies that the color-magnetic ob- 
jects should become lighter and more weakly coupled. 
They conjectured that electric components dominate the 
matter at high-T while the magnetic components domi- 
nate it near Tc- Furthermore, they gave arguments that 
an equilibrium point in which both electric and mag- 
netic coupling constants are equal (am = cte = 1) ex- 
ists at T lATc- Lattice observations have confirmed 
that monopoles form a relatively strongly coupled liquid 
where the magnetic coupling increases at high T (see 
discussion in J7]). Lattice observations of monopoles also 
strongly support the idea that they form a Bosc-Einstein 
condensate at exactly T ^ Tc ^TU\ . 

Recent two and three particle correlations in experi- 
ments at RHIC indicate that certain fluctuations occur- 
ring on top of the overall expanding quark-gluon plasma 
(QGP) have small or even zero expansion velocity in the 
near-Tc region, suggesting the existence of stabilized elec- 
tric flux tubes. In addition, the experimental results 
also show that the resulting clusters in AuAu collisions 
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are even larger than in pp collisions (i.e. in vacuum). 
The purpose of this paper is to use the classical dual- 
magnetohydrodynamic model to study their evolution 
when the dissipative effects are included. We will also 
attempt to evaluate the stability of the flux tubes in the 
near-Tc region, where we consider the plasma to be dual- 
magnetized 5 . Specifically, uncondensed color-magnetic 
monopoles circulate and form solenoidal currents which, 
by dual-Faraday's law, induce a color-electric field and 
create the flux tube. By studying the diffusion of this 
field in the medium we can calculate the halflifc of the 
flux tube and see if the field is strong enough to account 
for its stability during this magnetized phase. The dis- 
sipative effects are included by allowing the conductivity 
of the matter to be finite. 



B. Dual-Magnetohydrodynamics 

Magnetohydrodynamics (MHD) studies the dynamics 
of electrically conducting fluids and is widely used in 
plasma physics. In other words, it studies the interaction 
between a magnetic field and a plasma, treating it as a 
continuous medium. It is an approximation which keeps 
only the magnetic field in Maxwell's equations, while the 
electric field is entirely screened. Furthermore, MHD is 
used by solar physicists to describe the overall effects of 
electric currents and magnetic fields in the Sun's corona 
to describe the sun spots and flux tubes. 

As mentioned previously, one scenario has been pro- 
posed where the QCD matter near the region is dual- 
magnetized. In this region, the electric and magnetic 
screening masses are assumed to be Mm > Mg because 
the color-magnetic monopoles dominate the scene. (This 
is opposite to the case of high temperatures, in which 
electric particles - quarks and gluons - dominate, and 
therefore ^ g^T <C ^ gT, where g is the gauge 
coupling, small at high T. The high-T value for Mm kas 
been suggested by Polyakov [12]) It therefore is plausible 
to use DMHD ignoring the magnetic fields. 

In the usual plasmas like in the Sun, the electric screen- 
ing mass is very high and thus electric fields are ignored 
and magnetic fields are included. Under certain condi- 
tions the flux tubes can be formed as a result of solenoidal 
electric currents. Analogously, in RHIC collisions, there 
exist a "dual corona" ^ in which the electric flux tubes 
are produced by the magnetic currents. In what follows 
we study the diffusion of the field for flux tubes in MHD 
and then translate the results in the context of DMHD 
to investigate the lifetime of the flux tubes in the near-Tc 
region of the QCD matter. 



II. THE DIFFUSION EQUATION WITH FINITE 
CONDUCTIVITY 

The ideal MHD approximation is the limit of infinite 
conductivity of the plasma cr — > oo, analogous to the zero 



viscosity approximation in the case of ideal hydrodynam- 
ics. In this section, for self-consistency of the paper, we 
include the known derivation of the diffusion equation 
for the magnetic field in standard MHD with finite con- 
ductivity in order to account for dissipative effects. The 
fundamental equations of MHD are: 

f) Fi 

^ — =Vx£= (2.1) 



E^-vxB+- (2.2) 

and 

J=Vx^ (2.3) 

If the currents are orthogonal to the magnetic field (as 
they are in an ideal solenoid), then v x B = and the 
above equations simplify into the Diffusion Equation: 

§^,V'B (2.4, 

where we have used divB — and i] — I/Attct is the 
magnetic diffusivity. The magnetic diffusivity contains 
the conductivity of the plasma which will be discussed in 
part IV. The dual version of the diffusion equation is the 
same as the latter with the magnetic field B replaced by 
the dual-electric field E. In addition, instead of having 
ions and electrons with different masses (to^ >> m^), 
there are two (uncondensed) color-magnetic monopoles 
and antimonopoles with equal mass and opposite charge 
(7_i_ and g- forming two currents of equal magnitude cir- 
culating along the cylinder in opposite directions. 



III. SOLVING THE DIFFUSION EQUATION 

The flux tube is stabilized by magnetic currents circu- 
lating solenoidally on its surface. We place the solenoid 
in a cylindrical coordinate system with its vertical axis 
along the z direction. At t=0 each current ± J flows cir- 
cularly in the ±ip direction at r — Uq, the initial radius of 
the tube. The initial condition is taken to be a gaussian 
proflle of the field centered at the origin. At all times, the 
field points only along the z axis. We assume azimuthal 
symmetry so that there is no ip dependence and consider 
the field to be constant along the z direction. Therefore 
the field E{r,t) will only depend on the radial direction 
r and time t. The diffusion equation in cylindrical coor- 
dinates is: 

1 d ^^ dE{r,t) ^^ l dE{r,t) 
r dr dr rj dt 
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where 



V 



1 

Ana 



(3.2) 



is now the dual diffusivity and a is the dual conductivity 
of the plasma. The initial profile of the field is assumed 
to be a Gaussian centered at the origin: 



A{X) = XEo / e^^ Jo(Ar)rdr = - ^oAao^e^J ^'""^ 



Plugging this back in (3.4) we obtain an expression for 
E{r,t): 



E{r,0)=Eoe~ 



(3.3) 



Separation of variables E{r,t) — R{r)T{t) decouples the 
PDE (3.1) into two ordinary ODE's, one with time de- 
pendence and one with spatial dependence: 



dT{t) 
dt 



+ \^r^T{t) = 



The first equation gives us the time dependent part 



and the second equation is the Bessel equa- 



tion of zeroth order (m — 0) and the solution is a super- 
position of Bessel functions of the first and second kind. 
We exclude the Bessel function of second kind since they 
are not finite at r = and obtain R{t) — AJo(Ar). In 
Sturm-Liouville conditions the general solution is: 



E{r,t) 



AiX)Jo{Xr)e 



dX 



(3.4) 



To find the ^(A)'s lets look at E{r,t) at t = 0: 

E{r,Q) = / A(X)Jo(Xr)dX 
Jo 

Now multiply both sides by Jo{(3r)r and integrate over 
r: 



oo poo poo 

E{r,0)Jo{(3r)rdr ^ / A(X) / Jo{(3r)Jo{Xr)rdrdX 
^0 Jo 

The closure relation for Bessel functions gives: 



Jo(/3r)Jo(Ar)rdr = -5(/?- A) 
P 



and we obtain: 



r°° 1 1 
E{r,0)Jo{(3r)rdr = I A{X)-6{/3 - X)dX = -A{(3) 



A{X) ^ X / Eir,0)Jo{Xr)rdr 



E{r,t) 

and finally: 



\EaX flo^e" J Jo (Ar)rdr Jo(Ar)e-^'"*dA 
2 



Oo'' + ■irjt 



(3.5) 



Comparing (3.3) with (3.5), for t > Q, the term ac- 
quires an effective time dependent form + 4rjt that 
shows an increase in the field's squared radius. There- 
fore the field's time dependent radius is given by a{t) — 



IV. ESTIMATES OF THE CONDUCTIVITY 
AND THE FATE OF THE FLUX TUBES 

The standard electrodynamic plasma consists of elec- 
trons and ions having densities Ne and Ni respectively. 
Due to large difference in masses, mg ^ the current 
is assumed to be due to the motion of electrons only. 
The textbook expression for the conductivity of the usual 
plasmas, normal to the magnetic field [5^ is: 



0-_L 



for Z = 1 ions, and 



A-ne'^L^N, 

j^l/2y3/2 



(4.1) 



(4.2) 



which for E(r, 0) = Eo^- "° gives: 



is the collision rate between electron and ions. is the 
Coulomb logarithm and it is equal to ln(I/xmin) where 
Xrain is the magnitude of the smallest angles for which 
the scattering can still be regarded as Coulomb scatter- 
ing. The electron-electron collisions are ignored since 
they cannot change the total momentum of the electrons 
and thus do not modify the electron current. Our task is 
to translate these results into those appropriate for the 
QCD matter in the near-T'f. region. 

It consists of four main components. The first two 
are positively and negatively charged monopoles. Their 
densities are denoted by Nj^ and 7V_. As stressed before, 
the monopoles create counterdirected flows around the 
electric flux tube. The medium produced has overall zero 
magnetic charge, so that = iV_ = A'm/2. The last 
two components are electric objects - quarks and gluons 
- with densities and respectively. As noted in the 
Introduction, the electric components dominates at high 
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T, in the Quark-Gluon Plasma, but not in the near- Tc 
region, where they are suppressed. The total scattering 
rate for the positive monopoles is 



(4.3) 



and, by symmetry, the total collision rate is then vtot = 

It is important to note that the H — h scattering term 
is omitted for the same reason as the ee collisions are 
omitted in the electrodynamic plasmas: they would not 
change the total current. The H — cross section is the 
transport cross section due to magnetic Coulomb forces, 
so the H — collision rate is given by dual to the electron- 
ion rate above. It is obtained by the substitution of the 
coupling — (7^/47r (note the difference in 47r resulting 
from two different ways the fields are defined in QED and 
QCD) and the density N_. 

Scattering on electric objects is different, as discussed 
in details by Ratti and Shuryak [T3j. It has similar 
Rutherford-like scattering at small angles but the trans- 
port cross section is dominated by large (near-backward) 
scattering angles. In order to obtain expressions for the 
scattering rates and conductivity, we have to use certain 
empirical values of the parameters involved. All of them 
are, in principle, a function of the temperature T, but 
we restrict our discussion to the vicinity of T Tc- 
The values for the monopole density, magnetic Coulomb 
coupling and the mass that we use for the estimates are 
as follows: 



quantity 


value 


reference 




sa 1 


Fig.l of Rcf. MJ 




« 4/5 


Fig.3 of Ref.H 


m/Tc 


« 2 


Fig.6 of Ref.[IU] 




« 2 


Fig. 14 of Ref.'ni 



In addition, we assume that jz+g ~ v^q [22) . From 
these relations we obtain the -I— dual Coulomb collision 
rate: 



5, 



47rmi/2Ti'/2 
And from equation (4.3): 



47r2T? 



167r 



Vtot 



IGtt 



2Te + 2T, = T, 



giL,V2 + 647r 



The conductivity (4.1) is then: 



167r 



The value we used for Tc is the = 170 MeV from 
regular QCD with quarks rather than the larger = 260 
MeV value from pure gauge. Using (7^/47r w 4/5 and 



estimating 
get: 



= 1/10 so that L. 

167^LeT;^/2 _ 
25 ^ 

167rieV2 + 100 
25 



In(lO) w 2.30, we 
1113.04MeF = 5.65/m-^ 
) « 1793.04MeF = O.lO/m-^ 



3V27rTc 



25 ^ 
An I 



U.lSMeV = 6.16 x 10"^™"^ 



Finally, from equation (3.2): 



1 



6.56 X lO^^Mey^^ = 1.29/m 



Now we need to address the overall timing of the heavy 
ion collisions at RHIC. According to (very successful) 
hydrodynamical simulations, the duration of the magne- 
tized phase of the collisions is tm « 4 — 5 fm/c, at any 
centrality and any position of the fluid cell. According to 
our solution the mean square radius of the tube during 
this time is increasing by 



4?7rA/ 



(4.4) 



The energy per unit length e diffused during the magnetic 
phase is given by 



e{t) = 27r 



1 

8^ 



E{r,tf +K{r,t) 



rdr 



Where K — {l/2)pv'^ cx J{r,t)'^ is the kinetic energy 
density associated to the monopole current J(r, t). Using 
equation (2.3) we see that the time component of K scales 
like the time component of the magnetic energy. As a 
result, during the magnetic phase the total energy per 
unit length scales as : 



4:r]TM 



We estimate Oq « 0.5/m so that « 0.25/m^ and 
from the values above we get: 

23.26/to2 



M 



To gain more insight about the lifetime of the flux tube, 
let us calculate the half- life ti/2 of the field at the ori- 
gin. We know that at the origin and at i = we have 
£;(0, 0) = Eo- So, by definition: 



^(0,^1/2) 

and solving for ti^2' 



Eo 
2 



EoK 



ti/2 = T — '^'^at = 5 X 10 fm 
' 4 rj 



(4.5) 



(4.6) 



which is many times shorter than the expected lifetime 
of the observed flux tubes. 
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V. SUMMARY AND DISCUSSION 

In this paper we have used a classical dual- 
magnctohydrodynamic approach to calculate the flux 
tube lifetime in the magnetic phase of the QGP near 
the QCD phase transition. More specifically, we have 
found a solution for the flux tube including the dissi- 
pative "diffusive" term. We calculated the value for 
the "dual magnetic diffusion constant" using a picture 
of monopole-monopole and monopole-gluon rescattering 
and found that this crude classical rescattering model 
predicts very strong dissipative effects that are way too 
strong for the flux tubes to survive in the few fm/c time- 
frame of the magnetic phase. Yet, "ridge" correlations of 
the detected pions are found in experiments. 

One possible view on this, held e.g. by the BNL 
group [ini HI], is that the flux tubes do indeed decay 
very quickly, as the estimates above suggest, and the ob- 
served "ridge" is nothing but the extra amount of en- 
tropy left behind. A problem with this interpretation 
(discussed by one of us in [8 ) is that a spot of extra 
thermal entropy/energy would evolve hydrodynamically 
into a cylinder of several fm radius, which is in direct con- 
tradiction with the rather narrow (f) distribution width of 
the ridge. 

More likely, a classical approach to the flux tube dissi- 



pation is incorrect. First of all, unlike flux tubes usually 
considered by magnetohydrodynamics (e.g. in solar plas- 
mas) the QCD flux tubes under consideration are small in 
size and not larger than the quasiparticle Compton wave 
length. This suggests that one should use a quantum- 
mechanical description, like the one in [llj . Quantum 
effects in the monopole motion may provide two "super- 
currents" which propagate through each other without 
any dissipation. We know that this is the case in the 
confining vacuum. 

Finally, let us also mention that there are additional 
evidence for survival of the flux tubes in the magnetic 
phase, at T « Tc, which do not come from RHIC experi- 
ments but from lattice numerical simulations of the QCD 
thermodynamics. Those indicate that baryonic states re- 
main under such conditions, and that their density of 
state is well described by Hagedorn-like "stringy" states, 
see IHj. Unfortunately we do not know the lifetime of 
these baryons as the lattice thermodynamics does not 
give us such information. 
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